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Abstract 

We consider a dynamic market model where buyers and sellers submit limit orders. If at 
a given moment in time, the buyer is unable to complete his entire order due to the shortage 
of sell orders at the required limit price, the unmatched part of the order is recorded in the 
order book. Subsequently these buy unmatched orders may be matched with new incoming sell 
orders. The resulting demand curve constitutes the sole input to our model. The clearing price 
is then mechanically calculated using the market clearing condition. We use a Brownian sheet to 
model the demand curve, and provide some theoretical assumptions under which such a model 
is justified. 

Our main result is the proof that if there exists a unique equivalent martingale measure 
for the clearing price, then under some mild assumptions there is no arbitrage. We use the Ito- 
Wentzell formula to obtain that result, and also to characterize the dynamics of the demand curve 
and of the clearing price in the equivalent measure. We find that the volatility of the clearing 
price is (up to a stochastic factor) inversely proportional to the sum of buy and sell order flow 
density (evaluated at the clearing price), which confirms the intuition that volatility is inversely 
proportional to volume. We also demonstrate that our approach is implementable. We use real 
order book data and simulate option prices under a particularly simple parameterization of our 
model. 

The no-arbitrage conditions we obtain are applicable to a wide class of models, in the same 
way that the Heath- Jarrow-Morton conditions apply to a wide class of interest rate models. 



1 Introduction 

Most liquidity models in mathematical finance abstract the trading mechanism from the characteri- 
zation of prices in the resulting market. Our viewpoint is fundamentally different. In our model the 
equilibrium prices of the assets are completely determined by the order flow, which is viewed as an 
exogenous process. We model a market of assets without a specialist, where every trader submits 
limit orders, that is, for a buy order, the buyer specifies the maximum price, or the buy limit price, 
that he/she is willing to pay, and, for a sell order, the seller specifies the minimum price, or the sell 
limit price, at which he/she is willing to sellQ 

*With the collaboration of Thanh Hoang 
^Claremont McKenna College 
*Claremont Graduate University 

1 There is no loss of generality in this statement. A buy market order can be specified in our model as a buy limit 
order with the limit price equal to infinity. Since we model assets with only positive prices, a sell market order can 
be specified in our model as a sell limit order with a limit price equal to zero. 
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If at a given moment in time, the buyer is unable to complete his entire order due to the 
shortage of sell orders at the required limit price, the unmatched part of the order is recorded in 
the order book. A symmetric outcome follows in the case of incoming sell orders. Subsequently 
these buy unmatched orders may be matched with new incoming sell orders. We note that many 
electronic exchanges, such as NYSE Area [arcllj . operate like this. Time- priority is used to break 
indeterminacies of a match between an incoming buyer at at the limit price superior to the ask 
price, i.e., the lowest limit price in the sell order book. As the result, the equilibrium of clearing 
price process is always defined. 

Since the matching mechanism does not add any information to the economy, all information 
about asset prices is included in the order flow. Whether public exchanges should or should not 
reveal the real-time data contained in the order book is an important issue, which continues to 
preoccupy the financial markets community [WW02J. Our theoretical framework accommodates 
either viewpoint, but our empirical application is better tailored to the viewpoint that order books 
are public information. The current blossoming of the trading activity [EngOOj, [BLT06J, |H au08| . 
[AS08 , [BRZ09J seems to confirm our viewpoint that traders are (i) interested in understanding 
order book information, and (ii) trade on that information. 

We do not address in this paper the issue of differential information. The market microstructure 
literature (such as |Kyl85| , and all following models) considers various models of trading involving 
uninformed traders, called noise traders, and one or several informed traders. One of the key results 
of the Kyle model is that, given the information available to the noise traders, the resulting price 
process is a martingale in the appropriate measure, whereas it may not be for the informed traders. 
As a consequence we do not believe that abstracting issues of differential information is limiting. 
The order books reflect all the public information. We will show (under certain conditions) that 
the clearing price is a martingale in the filtration corresponding to the public information, and not 
private information. 

There are roughly two different class of models in the liquidity literature. The first class of 
models ( |.Tar92j . [Jar 94] . |PS98bj . |PS98aj . |Fre98j . |SW00| . [BB04| . [RSTD] ) considers the action of 
a large trader who can manipulate the prices in the market. There are mainly two different types of 
strategies a large trader can employ to that effect. The first one is to corner the market, and then 
squeeze the shorts. The second one is to "front-run one's own trades". While some exchanges have 
rules to curtail the cornering of the market, the front-running seems more difficult to ban from an 
exchange. It is known in discrete-time trading, that, if there is no possibility of arbitrage for small 
traders in periods where the large investor does not trade, then there is no market manipulation 
strategy. In this paper, we assume that Jarrow's conditions ( | Jar94] ) for the absence of the market 
manipulation strategy in discrete time hold. We exploit then the theoretical results from [BB04 , 
and |KR09| to prove the absence of arbitrage in our continuous-time model, under certain conditions. 
The second class of models ( |C/ JP04| , |QR07| , [QST10|, |GSllj abstracts the issues of the market 



manipulation away, and considers all traders as price-takers. In particular, QJP04 introduced an 
exogenous residual supply curve against which an investor trades. The investor trades market 
orders, and his/her order is matched instantaneously. As a consequence of the instantaneity, it is 
plausible for [QJP04| to assume the "price effect of an order is limited to the very moment when 
the order is placed in the market" (dixit [BB04]), so that that the residual supply curve at a future 
time is statistically independent from the order just matched. For us however, since all information 
is contained in the order flow, this assumption is not plausible, since it does not explain how prices 
can incorporate information from the arrival of new orders. 

The key to our approach is to distinguish between what we call the cross orders, i.e., the orders 
submitted at a limit price such that they are likely to be instantaneously matched, and what we call 
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the uncross orders, i.e., orders which will spend a positive amount of time in the order books before 
being either matched or cancelled. We are not aware of workable assumptions in the literature that 
warrant the absence of arbitrage strategy involving the uncross orders. For cross orders however 
the task is simpler, and our approach takes advantage of the results presented in |BB04j . which 
assume that trading is immediate. Anecdotal evidence shows however that market manipulation 
strategies tend in practice to occur over short periods of time, casting doubts on the practicality of 
implementing market manipulation strategies with uncross orders. 

In our model all the information is contained in a Brownian sheet. The Brownian sheet drives 
the dynamics of the demand curve. The advantage of using a Brownian sheet is that we have the 
same cardinality of independent sources of noise (namely, the cardinality of a real interval) as the 
cardinality of the set of exogenous stochastic processes. Our methodology generalizes the nonlinear 
partial differential equation approaches contained in the literature cited above, in the same way 
that the Heath-Jarrow-Morton methodology generalizes term structure models. 

The main result in our article is to prove that if there exists a unique equivalent martingale 
measure Q for the clearing price, then (under some mild conditions) there is no arbitrage. We 
use the Ito-Wentzell formula to obtain that result, and also to characterize the dynamics of the 
demand curve and of the clearing price in the measure Q. We find that the volatility of the clearing 
price is (up to a stochastic factor) inversely proportional to the sum of the buy and sell order 
flow densities (evaluated at the clearing price), which confirms the intuition that the volatility is 
inversely proportional to the volume. We also demonstrate that our approach is implementable. 
Although we do not prove a second fundamental theorem of asset pricing, we naively use a special 
parameterization of our model to price options. As in the early days of the Heath-Jarrow-Morton 
methodology, we use only historical estimation (in our case, of the order book) to fit our model and 
to solve for the market price of risk. We expect that, should this paper meet with interest around 
practitioners, market-implied implementations will see the day. Unsurprisingly, we obtain a smile 
curve for the implied volatility. We note that this particular feature is not a very strong sign of the 
adequacy of our approach to model asset prices, as most models that came after |BS73j result in a 
smile curve for the implied volatility. Although limited, our results are however encouraging. They 
show that a fairly demanding theoretical model can be easily implemented. 

The structure of the paper is as follows. In Section [2j we introduce two classes of models, one 
with atomistic traders only, and one with atomistic traders and a large trader. In the first one, 
almost by definition, the demand curve turns out to be continuous in time. In the second one, 
continuity (under a set of assumptions) is proved. For both models, we prove that there exist one 
or several martingale measures Q for the clearing price and that under these measures there is no 
arbitrage. In Section [3j we characterize the price process under the risk- neutral measure, giving the 
conditions under which it is unique. In Section [4j we describe our data set and the methodology we 
used to extract the relevant information. We also present our implementation, namely a simulation 
of option prices under the risk- neutral measure. 

2 Model 

We assume a filtered probability space (f2, T, {Ft}, IP). We will see later what processes generate 
the filtration. 
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2.1 The Market Mechanism 



A buy limit order specifies how many shares a trader wants to buy, and at what maximum price 
he is willing to buy them. We call this price the (buy) limit price. A buy limit order specifies how 
many shares a trader wants to buy, and at what maximum price he is willing to buy them. We call 
this price the (sell) limit price. Both buy and sell limit prices are denoted by p, and should not be 
confused with the clearing price, which is denoted by 7r(t). The unmatched buy and sell orders are 
kept in order books until they are either cancelled or matched with an incoming order. An incoming 
order is matched with the order in the opposite side of the market which has the best price. The 
clearing price of the transaction is equal to the limit price of the order in the book, and not of the 
incoming order. Partial execution is allowed, and ties are resolved by the time-priority. Below we 
present an example of the matching mechanism in discrete time, that is at most one order arrives 
at time t £ {0,1,2,..}. 

Example 2.1. Suppose that the clearing price at time is any price 7r(0) G [100,120]. After 
clearing, that is, when < t < 1 we suppose that the order book contains the following orders 



Buy Order Book 
Price I Quantity 



100 



10 



Sell O 
Price 


rder Book 
Quantity 


120 


10 


130 


10 



At time t = 1 a buy order arrives with a limit price of $125, and a quantity of 15. The exchange 
matches it with the best sell order, i.e., the one with a sell limit price of $120. However, the 
execution is only partial, and the remainder of the buy order is placed in the order book at the limit 
price of $120, resulting in the following order book: 



Buy 
Price 


rder Book 
Quantity 


100 


10 


125 


5 



Sell Order Book 


Price 


Quantity 


130 


10 



The clearing price at time 1 is equal to the limit price of the sell order, i.e.: 



tt(1) = 120. 

This example illustrates several properties of the limit order markets. First, the clearing price 
is always defined, and can assume any positive value [^J 

Second, it is not inconceivable that an incoming order "crosses" the order book, i.e., for the case 
of a buy order, that it best higher limit price than the best sell order limit price, or best ask, since 
the buyer does not lose a cent. Crossing the book is indeed advantageous for two reasons: first, 
it allows for faster execution. In our example, had the buyer submitted an order at price $130 he 
would have bought the complete quantity of shares (15) that he desired, rather than waiting an 
indeterminate amount of time until enough sell orders arrive at his limit price. Second, suppose 
that several buy orders are submitted at the same time. In case the demand exceeds the supply at 
the best ask, the buy orders with the highest limit price are executed first. Our own data analysis 
(see section [4]) shows that few orders cross the NYSE ArcaBook [arcllj . This is consistent with the 
theory of optimal order book placement suggested by Rosu |Ros09j . 

2 We do not consider markets for swaps, where the prices can be negative. 
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Assumption 2.1. Buy and sell limit prices can assume any real value between and S. They are 
usually denoted by p. Orders can be submitted to the market at any time t £ R + . 



2.2 The Brownian Sheet 

We now turn our attention to the continuous-time case. Note that we do not prove convergence 
of a discrete-time model to our continuous-time model. We take the latter as a given, plausible 
model of the market. We start with a probability space (Q,J-,¥) satisfying the usual conditions. 
The uncertainty is described by a one-dimensional Brownian sheet W(t, s), which is a continuous 
version of the construction (4.20) in [DPZ92J p. 100, i.e.: 



Here {(3j} is a family of independent real- valued standard Wiener processes, and {gj} is an 
orthonormal and complete basis for some Hilbert space. The variable t refers to time and the 
variable s identifies the "factor" information necessary to model a large collection of processes 
identified by the variable p, with < p < S. While the bound S is usually taken to be equal to 
one in the literature, we assume instead that S is a large value for ease of modeling, as will become 
clear later. The filtration {J-t} is generated by the collection of Wiener processes {/3j}. A stochastic 
integral of an JVadapted integrand o~(t, s) with respect to the Brownian sheet is denoted by: 



2.2.1 A Market with Atomistic Traders 

Definition 2.1. The net demand curve Q is a function [0, S] x M + x fi— > R, which value Q(p, t, u) 
is equal to the difference between the quantity of shares available for purchase and the quantity 
of shares available for sale at price p and at time t. For each p the stochastic process Q(.,t, .) is 
Ft-adapted. The random variables Q(p,t, .) are assumed to be uniformly bounded in p and t almost 
surely. 

The net demand curve represents then the information available in the (limit) order books of 
the exchange. 

Remark 2.1. As we will see, the process Q will be assumed to be continuous in time, and we will 
prove that the clearing price is continuous in this model. Thus we do not need to specify in this 
model whether Q represents the net demand just before clearing, or right after clearing, since these 
quantities move continuously. The same remark will not apply a priori in our model with a large 
trader. 

Remark 2.2. The net demand curve is decreasing in p. 

Definition 2.2. The clearing price ir(t) is an J-t-adapted stochastic process which either satisfies 



when there is a solution to Q), or is otherwise defined by continuation, i.e., n(t) is equal to the 
value of ir at the latest time s <t for which there was a solution to Q(ir(s-), s) = 0. 




0<a<s 





(1) 
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When the net demand curve is continuous in time, we have of course, for any t: 

Q(ir(t),t) = Q(w(t),t+). (2) 
Definition 2.3. A limit order submitted at time t crosses the market at time t if: 

• either it is a buy order with limit price p > 7r(i_) 

• or it is a sell order with limit price p < 7r(t_) 

We call these orders cross orders. All the other orders are called uncross orders. 

In practice it is rare that limit orders cross the market, indeed there are equilibrium models in 
which such orders are not rational [Ros09j. 

Remark 2.3. There is a more general way to model a market. We could accumulate all of the buy 
order quantities with limit price higher than p that enter the system (and withdraw the cancelled 
quantities when an order is cancelled) into a cumulative demand curve T>(p,t). Likewise we could 
accumulate all of the sell order quantities with limit price less than p that enter the system (and 
also withdraw the cancelled quantities) into a cumulative supply curve S(p, t). We call these curves 
"cumulative" because we include the order quantities that are matched in them, and thus, if no 
order is cancelled, these curves increase in t. We observe that: 

Q = V-S 

Thus, modeling only the net demand curve is reductive. While there is a (not very natural) way 
to extend the results of this paper to a model including both the cumulated demand and supply 
curves, we decided to present here only our simpler model for the following reasons. As we argued 
in the previous section, orders rarely "cross", there is a numerical instability when modeling both 
T> and S. 

Remark 2.4. Due to Assumption |2.4[ the clearing price is limited to take values between and 
S. The frontier ir(t) = corresponds to a bankruptcy, and the frontier ir(t) = S corresponds to a 
higher limit (say S = $1M) set by the exchange to prevent excessive speculation. 

Remark 2.5. Questions of uniqueness of the clearing price will be addressed when we describe the 
stochastic differential equation that the clearing price (should it exist) satisfies. 

The next assumption is standard. Note that by "transaction costs" we do not mean the liquidity 
cost incurred, but an additional cost per transaction that the exchange would charge the traders. 
For instance, a buyer would pay per share a cost 7r(t) + c, with c > 0. 

Assumption 2.2. The market is frictionless, i.e. c = 0. 

We now develop a model for the order quantity specified in any order. Modeling a net demand 
curve which is twice differentiable with respect to price is clearly easier than modeling a discrete 
curve, and we assume it. It can occur when there is an uncountably infinite amount of traders, and 
traders are atomistic. We need however to rule out a degenerate case when a subset of non-zero 
measure of all the traders agree on the limit price, thus generating a discontinuity in the demand. 
This can be justified by assuming differential information amongst traders about the real value of 
the stock, as in the market microstructure literature. Note that the atomistic traders do not need 
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to be noise traders for the market to be consistent. Indeed, we can (but do not need to) assume that 
all traders know all the information contained in the order books at all times. If they do know this 
information, as specified in the introduction, the clearing price does not bring any extra information 
as the traders can compute the clearing price by themselves at all times [^J Making the order book 
public information can be an important assumption in some settings. Indeed, our empirical studies 
(see Section [4]) show that the order flow is tightly concentrated around the current clearing price. 
This confirms the obvious intuition that traders are strategic and choose their limit price based 
on their (slightly delayed) knowledge of the clearing price. In this paper we will fit a parametric 
model to the market which exploits this intuition. Our model is fairly robust when the net demand 
is modelled as a function of the difference in price p — 7r(i), but not robust when it is modelled 
(as discussed above) as a function of the limit price p without a reference to the whole order book 
information. 

Assumption 2.3. There is a continuum of atomistic buyers and sellers who trade on the market. 
The resulting net demand curve Q is twice differentiable in price p and continuous in t. We assume 
that 



dQ 
dp 
dQ 
dp 



_ dQ 

< 0, for < p < S 



= 

p=S 



p 

For the moment we just state the stochastic differential equation that Q satisfies, assuming that 
the resulting net demand curve satisfies Assumption |2,3[ We defer the task of showing examples 
where these conditions are satisfied to the next chapter, since in the next model the exact same 
conditions will have to apply, and it is more convenient for the reader to have a single place with 
the exact specification of the model. Thus 

dQ(p,t) = n Q (p,t)dt + o- Q (p,t) [ b Q (p,s,t)W(ds,dt), for 0<p<S (3) 

Js=0 

Q(P,0) = Q (p), for 0<p<S 

The coefficients pq, o~q, and 6q are J-t-adapted. For the moment we just assume that they are 
such that the solution to ^ exists, is unique, and is uniformly bounded in p and t almost surely. 
Besides, for every p, the process Q(p, ■) is a semimartingale. Also, we enforce: 

S 



/bq(p, s,t)ds = 1, for every p and t. 



Definition 2.4. A (trading) strategy = (0(t)) is a semimartingale that represents a number of 
shares held by the investor at each point in time. If the strategy is self -financing (see e.g. [BBOJ^] 
for a definition) the process f3 e representing the value of the cash account is uniquely defined. 

We refer the reader to |Jar94j for a definition of market manipulation strategies in discrete time. 

Definition 2.5. For every real-valued x an inverse process P(x,t) satisfies 

Q(P(x,t),t)=x. (4) 

The process P(x, t) is undefined when Q) does not admit a solution. 

3 In order to account for the delay in information processing, it would be more plausible to assume that the 
information available to all traders at time t corresponds to Ft- and not Tt but this would complicate the model 
significantly, while yielding not much extra conceptual value. 
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Remark 2.6. Since Q is strictly monotonic in p, then whenever Q exists, it is also unique. 
Definition 2.6 ((3.1) in [BB04]). The asymptotic liquidation proceeds L($,t) are defined as: 



L(&,t) = [ P(x,t)dx. 
Jo 



This definition is from [BB04] for the proceeds of a fast liquidation strategy of a large trader 
from i? to 0. The intuition behind this process will become more clear when we consider a market 
with a large trader. 

Definition 2.7. The real wealth process achieved by a self- financing trading strategy 9 is given by 

V e (t) = (3 e (t) + L(9(t),t) 

Lemma 2.1 (Lemma 3.2 in [BB04]). For any self-financing semimartingale strategy 9, the dynamics 
of the real wealth process V 8 are given by 

V e {t)-V e (0-) (5) 

= f L{9{u-),du)-\ f P'(9(u^),u)d[9,9} c s - V f ^ {P(9(u),u) - P(x,u)}dx. 
Jo ^Jo o< v <tM»-) 

Definition 2.8. An arbitrage (strategy) is a self-financing trading strategy 9 such that V s (0-) = 
and 

F{V e (t) > 0) > 0, 
F(V e {t) > 0) > 0. 

Theorem 2.1. Suppose in addition to our standing assumptions that 

CI) for self -financing strategies involving only cross orders, J arrow's jJar94f discrete-time con- 
ditions for absence of market manipulation strategy hold, 

C2) no arbitrage strategy involves uncross orders, 

C3) the volatility o~Q(p,t) is bounded away from zero, uniformly in p, 

C4) there is no path such that Q(S,t) > or Q(0,t) < 0. 
Then 

Fl) there exists at least one martingale measure Q for J L($,dt), 

F2) there is no arbitrage strategy, 

F3) the clearing price n(t) is continuous, 

F4) any such measure Q is also a martingale measure for n(t). 



Proof. All arguments of the proof of Theorem 2.2 can be used, replacing the (non-infinitesimal) 
strategy 9 of the large trader by an infinitesimal strategy 9dp of an atomistic trader. Since Theorem 



2.2 is more general, we will only show the proof of Theorem |2.2[ □ 



A market composed only of atomistic traders is not very realistic. In the next section, we will 
show that under certain conditions a large trader can be added to our market, and the no-arbitrage 
conditions can still be obtained. 
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2.3 A Market with Atomistic Traders and a Large Trader 

A continuous (in time) net supply curve can also arise when large traders decide to split their 
large orders into atomistic orders. CJP04 show that this is indeed an optimal strategy in their 
model. The same decision turns out to be optimal in our model, under certain conditions which 



are elaborated in the next assumptions. We will prove this fact in Theorem 2.2 thus making our 
model is self-consistent. 

Compared to |CJP04| , there are two main conceptual difficulties when adding a large trader 
into the model. First, the trader may submit orders which are not instantaneously matched, unlike 
the market orders in the |CJP04| model. This is why we consider separately in this section the 
cross orders and the uncross orders. The cross orders will be matched instantaneously, and thus 
they behave like market orders in the [CJP04| model, and continuity (in time) can be proved. For 
uncross orders, one would need a more detailed economic model to specify in which cases it is 
advantageous for traders to submit a continuous (in time) order flow. Rather than going into these 
details, we assume that this holds true. In any event, market manipulation strategies (such as 
market cornering) are probably more likely to be implemented with cross orders than with uncross 
orders. 

The second difficulty is that we have to prove that the large trader cannot manipulate the market. 
The conditions in [Jar 94] apply when traders can submit orders at discrete time intervals. [BB04J 
and |KR09j show the conditions under which the discrete-time conditions extend to continuous time, 
and our proof of consistency of the market with a large trader will consist in part in checking that 
[KR09J conditions apply. 



Since the continuity in time of the net demand is not assumed any more (see Remark 2.1), we 
need now to distinguish between the incoming orders (or the order flow) of the large trader, and 
the order book position. 

Definition 2.9. The net demand curves of a large (atomistic) trader Ql (Qa) is a function [0,P] x 
]R + x S7— > R whose value Qi(p,t,uj) (QA(p,t,co)) is equal to the difference between the quantity of 
shares submitted for purchase and the quantity of shares submitted for sale at price p at time 
t. For each p the stochastic processes Qz,(.,t, , ) and QA(-,t, ■) are Tt-adapted semimartingales. As 
before the net demand of the atomistic traders satisfies 

dQA(p,t) = fj,Q A (p,t)dt + o-Q A (p,t) f b QA (p,s,t)W(ds,dt) forO<p<S, (6) 

Js=0 

Qa(p,0) = Qa,o(p) for 0<p<S. 

Definition 2.10. For every real-valued x the process Pa{x,i) satisfies 

Q A {PA(x,t),t)=x. (7) 

Definition 2.11. The asymptotic liquidation proceeds of the large trader Li,($,t) are defined by 

ri9 



rv 

L L (0,t) = / P A (x,t)dx. 
J 



Assumption 2.4. Both Ql and Qa are twice differ entiable in p. Only Qa is assumed to be 
continuous in t. 

Remark 2.7. The (total) net demand curve satisfies 

Q = Ql + Qa- 
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Assumption 2.5. For simplicity we assume 



Q(0,t) > 0. 

Assumption 2.6. For each p > ir(t), the function i) is continuous in time. 

Theorem 2.2. Suppose in addition to the standing assumptions that 

CI) for self-financing strategies involving only cross orders, Jarrow's ^JarU^j discrete time con- 
ditions for absence of market manipulation strategy hold, 

C2) no arbitrage strategy involves uncross orders, 

C3) the volatility aQ A (p,t) is bounded away from zero, uniformly in p, 
C4) there is no path such that Q(S,t) > or Q(0,t) < 0. 
Then 

Fl) there exists at least one martingale measure Q for j Ll($> dt), 

F2) there is no arbitrage strategy, 

F3) the net demand curve Q is continuous in t, 

F4) the clearing price 7r(t) is continuous, 

F5) any such measure Q is also a martingale measure for ir(t). 
The full proof is presented in the Appendix. Here we show only a summary of the proof. 

Summary of the Proof: We verify that the |KR09] conditions hold. Thus j L^-Q^t) is a 
Q-martingale, and no market manipulation exists that involves cross orders. We prove that tame 
strategies (9 continuous in t) are optimal for the large trader so that, first Q is continuous in t., 
and, second, cross orders are traded only at the clearing price ir(t). 

□ 

Remark 2.8. Under certain conditions it is possible to extend our model to more than one large 
trader on the market. Following the assumptions of Theorem |2.2[ the net supply curve Q in a 
market with one large trader is indistinguishable from the net supply curve in a market with only 
atomistic traders. It is thus plausible that, should a second large trader arrive in the market, he 
would behave like the first large trader. However, there are many different ways to justify this 
result, and we believe that this discussion would be more appropriate in an economics journal than 
here. 

3 Characterization of the Price Process in the Risk-Neutral Mea- 
sure 

In this section we characterize the price process in either one of the market models we specified 
earlier, since they turned out to be equivalent for our purposes. 

Assumption 3.1. To avoid repetition we assume in this section that there is no path such that 
Q(S,t) > or Q(0,t) < 0. 

Definition 3.1. The market price of risk X is a function [0, S] x M + x £1— > ]R. The market price of 
risk process A(s, ., .) is an Tt-adapted semimartingale for every s £ [0, S]. We define the <Q-measure 
as a measure such that the process is a Brownian sheet, where: 

W q (ds,dt) = W{ds,dt) + X(s,t)dt. (8) 
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We can then define the clearing price process as 



dir(t) = a w (t) / b 7T (s,t)W Q (ds,dt), 



with 



b^{s 1 t) 2 ds = 1. 



(9) 



Since Q(p,t) must be strictly decreasing in p, we find it convenient to slightly modify the 
definition into: 

Q(p,t) = Q(p,0)- [ P q(y,t)dy, 
Jo 

where we define Q(0, t) and q(p, t) (for < p < S) to be strictly positive processes with 

dQ(p,t) = HQ(0,t)dt-<r Q (0,t) jb q (0,s,t)W(ds,dt) Q(0, 0) = Q o (0), (10) 
dq(p,t) = n q (p,t)dt + a q (p,t) b q (p,s,t)W{ds,dt) q(p, 0) = Q (p), (11) 

is 

for < p < S 
q(0,t) = 0. 



Like before, the coefficients of the equations (10) and (11) are J^-adapted. We assume that the 



solution to (|10|) and (| 1 1 p exists, is unique, and is uniformly bounded in p and t almost surely. Also, 

"5 



we require that 



/ b^(p, s, t)ds = 1 for every p and t. 
Js=0 

Remark 3.1. The process q is a density of orders. By definition 

q(p)dp = quantity of shares available for purchase with limit price in \p,p + dp] 
+ quantity of shares available for sale with limit price in [p, p + dp] . 

Remark 3.2. Assuming that Q(0,t) is twice-differentiable in p, we must make sure that q(p,t) is 
differentiable in p, for the process Q(p, t) to be twice-differentiable in p. This occurs if, for instance, 



dq(p, t) 



(p- s)W(ds,dt). 



s=0 



We now define the following processes: 

C(ir,t) = -a n (t) (J?- (a Q {0,t) ^6 ? (0,s,i)6 7r (s,t)dsj + <r q (ir(t), t) J b,,( tt. s. I )/,_(*. / ),/. s 

f n 1 da 

b(7T,t) = -p Q (0,t) + J fl q ( P ,t)dpdt+--^(7T,t)(a 7I (t)) 2 -C(7T,t), 



f"K 

JO 



(p, s, t)ds 
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Remark 3.3. If one is not interested in modeling the correlation between orders at a different 
limit price, one may think that a Brownian sheet is not necessary, in the same way that the Heath- 
Jarrow-Morton model is often implemented with only 2 factor^} As we shall see however, the full 
complexity of a Brownian sheet is necessary for the market price of risk equations to have a solution. 

Definition 3.2. The market price of risk equations are: 




E(ir,s,t)X(s,t)ds = b(TT,t), for < vr < P. 



Theorem 3.1. Suppose that the previous assumptions hold true. In addition, suppose that the 
market price of risk equations have a unique solution. Then there is no arbitrage. 

Proof. A market clears if Q(p(t),t) = 0, or, equivalently, if dQ(p(t),t) = 0. We use the Ito-Wentzell 
formula to compute dQ(p(t),t) and set dQ(p(t),t) = 0. 



fi Q (0,t)dt 



0+ 



/i q (p, t)dpdt 



cr q (p,t) / b q (p,s,t)W(ds,dt)dp 



(12) 



q(ir(t),t)a n (t) / b n (s,t)W (ds,dt) 



ldq 
2dp 



Tr(t),t)(a n (t)) 2 dt + C(Tr(t),t)dt = 0. 



We equate the volatility terms to zero above and find 

Jo Vq(p,t)b q {p,s,t)dp 



a n (t)b n (s,t) 



q«t),t) 



(13) 



We substitute (|8j) in (12) and equate the drift terms to zero. This results in 

rw(t) 



s JO 



o~ q (p, t)b q (p, s, t)dpX(s, t)ds 

,7T(*) 



ldq 



fiq(p,t)dpdt+ ^(vr(t),t)K(t)) 2 - C(n(t),t). 



-mq(o,*) + 

Since the above must hold for any value of vr(t), then the market price equations must be satisfied. 
Thus there exists a unique measure Q such that n is a martingale. By theorems 2.1 and 2.2, there 
exists a non-empty set Q of martingale measures for J L($, dt). Besides, any measure Q G Q must 
be a martingale measure for n. Uniqueness of Q ensures that Q = {Q}, thus Q is a martingale 
measure for J L($,dt). Therefore Theorems 



2.1 



and 2.2 imply no arbitrage. 



□ 



Remark 3.4. In a numerical implementation, the market price of risk equations will be a set of S 
linear equations with S linear unknowns. Generically, like other market price of risk equations in 
finance, these will admit a unique solution. This was the case in all our simulations. 

Remark 3.5. Integrating (13) and using Q, we see that 

(Io S (Io {t) ° q {p,t)b q (p,s,t)fdpds)V 2 



a n (t) 



(n(t),t) 



(14) 



The denominator of ( 14 ) shows that more orders at the clearing price decrease volatility. This is to 
be expected. The effect of the numerator is harder to analyze, and shows that the volatility of the 
whole net demand curve affects the volatility of the clearing price. 



4 Note however that Carmona and Tehranchi [CT06| show that the Heath- Jarrow-Morton model, if there are less 
factors than forward rates, may result in a special type of arbitrage. 
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Before we move to the implementation of our model, we perform an empirical analysis of the 
market, which will guide us in specifying a parametric model. 



4 Empirical Analysis 

4.1 Implementation 

To the best of our knowledge, our methodology to model prices is quite original. In the same way that 
the HJM methodology opened the way to the development of several parametric models, we hope 
that this paper will result in an effort to develop models for the demand curve that are appropriate 
for risk management. This section only scratches the surface of that effort. Alternatively, this is 
what practitioners would call a "proof of concept", showing that the model can be implemented 
without too much effort. We do not make any claims about the power of the model we present 
hereafter, which is perhaps the simplest no-arbitrage model imaginable in our framework. 

Compared to the Heath-Jarrow-Morton model for the forward rates, it is necessary in our model 
to know the drift of Q(p, t) in the physical measure in order to simulate the price in the risk- neutral 
measure. In practice, there are two different methods to determine this drift: historical estimation 
or market implied. In the historical estimation method, one calibrates the physical drift /iq(0, t) in 
the model to market observables, and then solves for the market price of risk. In the market implied 
method, one calibrates the model to, say, a smile curve of option prices, like for other stochastic 
volatility models. This paper fell short of proving a second fundamental theorem of asset pricing, 
which would justify the implied market method. This will be a goal a subsequent paper. 

In our implementation, we first use the historical estimation method to estimate the drift of Q 
in the physical measure, and then calculate options prices by simulation, in essence assuming that 
the second fundamental theorem works. 

4.2 A Model with Relative Prices 

As explained in the introduction, and shown later, a model will be more robust if it assumes a 
direct dependence between the net demand curve Q(p, t) and the relative prices p — ir(t). This takes 
into account the fact that the investors observe the market in real-time. We are then obligated to 
introduce Q(p, it, t) as the net demand curve at price p when the clearing price is tt. We have then 
as before 



r (t) 



Q(p,7T(t),t) 



Q(0,vr(t),t)- / q(p,n(t),t)dp, 



(15) 



dQ{0,Tr(t),t) 



dq{p,ir(t),t) 



H Q {Q,Tr{t),t)dt - <7,(0,7r(t),t) / b,(0, vr(t), s, t)W(ds, dt). 



H q (p,ir(t),t)dt + a q (p,ir(t),t) / b q (p, ir(t), s, t)W(ds, dt), 




Let 




o 



Then 



dQ{p,TT,t) 

dir 



- q 



(vr,7r,t) - / —q(p,7r,t)dp 




o 



13 



and 



„ , , d 2 Q(0,ir,t) . . . 8 . f d 2 . 

F^{tt, t) = q p {TT, tt, t) - q % {ir, tt, t) - T^gi^, x,t) ~ / 3^2 9 7r ' f ) d P- 



We define 



TT 

H(s, 7T, t) = ~ J a q (p, TT, t)b q (p, TT, s, t)dp. 



TT[t) 

dQ{TT{t),t) = fi Q (0,7T(t),t)- f ^ q { P ,TT(t),t)dt-O q { Pl TT{t),t)(fb q {p,TT(t),S,t)W{ds,dt) 



p=0 

Hence 

d f Q(7 

— H(s,TT,t) = - ^[a q (p,TT,t)b q (p,TT(t),S,t)}dp- a q (TT,TT,t)b q (TT,TT,S,t). 
p=0 

Therefore the Ito-Wentzell's formula takes the form 

dp 

\J s / 

p=0 

+F n (TT(t),t)dTT(t) + -F^(7t(t), t)(dTT) 2 + C(n(t),t)dt, 

where 

— H( S ,TT(t),t)ds. 

Market clears if dQ(ir(t),t) = 0. Equating the volatilities results in 

f Jo (7 q (p,TT{t),t)b q {p,TT(t),S,t)dp 

= E-^J) • 

Apart from the denominator, this is the same equation as ( |13| ). We now define 

b{TT,t) = J N {p,TT,t)dp- ll Q {Q,TT,t) - ^F n7T (TT(t),t)al(t) - C(ic(t),t). (16) 

Finally, the market price of risk equations can be expressed as 

4.3 A Parametric Implementation with Relative Prices 

We discretize our model in 3 dimensions: time t, limit price p and factor s. For the market price 
of risk equations to have a unique solution the number of price buckets where orders are assigned 
should be equal to S, namely the number of factors. We let Ap be the size of a price bucket. 
We define the relative price k as 

k = p — TT, 

and the relative net demand Q(k,TT,t) as 

Q(k,TT,t) = Q{p,TT,t). 
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We simplify the model above (as expressed in (15)), and assume that the order flow depends 
only on the relative price k, i.e., Q(k, tt, t) = Q(k, t). Now assume without loss of generality that S 
is even, and define K = 5/2. We model the relative net demand curve as 



k 

Q(k,t)=Q(-K,t)+ £ ?(ht). 

l=-K+l 

In our model, the logarithm of the order flow quantities follow Ornstein-Uhlenbeck processes. 
This ensures the stationarity as well as the positivity of the order flow quantities. In other terms 

K 

dlog Q(0 ) t) = -o (0)(logQ(0,t)- log Q(0))dt + a r Q el (0) ]T b q (k,j)^/ApdW 3 (t) (18) 

j=-K+l 
K 

dlogq(k,t) = -a q (k) (log q(k,t) -logq(k))dt + a r q el (k) ^ b q (k,j) y /Ap~dW j (t), with (19) 

j=-K+l 

k = -K + 1..K 

It is thus necessary to estimate the parameters a q (k) > 0, ag(0) > 0, (T q el (k) and b q (k,j), as 
well as the initial values. We notice that our model is not twice differentiable in p, however, for 
our parameter values, it does not behave significantly differently from a smoothed version of that 
model. 

4.3.1 Results 

Data: We collected high frequency data from NYSE ArcaBook |arcllj for General Electric Com- 
pany (GE) on April 1 st , 2011. GE limit order is characterized by six intrinsic quantities: (1) message 
type: e.g. "A": add new order, "M": modify order, "D": delete order; (2) trading type: e.g. "B": 
buy limit order or "S": sell limit order; (3) time: recorded when a new event occurs, e.g. add 
new order, modify order or delete order; (4) ID: a unique identifier of each limit order; (5) price in 
dollars; and (6) size in number of shares. 

We selected data during the trading time from 9:30 AM to 4:00 PM EST. Before doing any 
analysis, we removed obvious errors in the data, e.g. abnormal prices of USD 122 or USD 0.01. 
We managed to keep track of approximately 90% of the original limit orders. For the cancellation 
of limit orders, we assumed that if the amount of time was less than or equal 2 minutes when the 
"delete" message occurred after the "modify" message of the same order, this order was cancelled. 

Let p m in be the minimum price of USD 20.00 and p max be the maximum price of USD 20.62 
during the trading day. We partitioned the price space into 5 cents bins 

p min = 20.00 = P-K < P-K+l < ■ ■ < PK-l <PK = Pmax = 20.62 

PK+i-PK = Ap (5 cent). 
Likewise, we partitioned the time into: 

tj+i — ti = At (1 minute) 

Clearing Prices ir(t): By dividing the time period from 9:30 AM to 4:00 PM EST into one- minute 
time intervals, we obtained 390 one- minute time intervals (6.5 trading hours * 60 minutes). 
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Table 1: Statistical Summary for ir(t) 



Summary 


Output 


Summary 


Output 


nobs 


Ton 


oil/ iviean 


n nn/i 1 7 a 


NAs 


o 


LCL Mean 


20.3294 


Minimum 


20.05 


UCL Mean 


20.3458 


Maximum 


20.61 


Variance 


0.006795 


1. Quartile 


20.29 


Stdev 


0.082435 


3. Quartile 


20.39 


Skewness 


-0.289841 


Mean 


20.338 


Kurtosis 


0.277282 


Median 


20.340 


Sum 


7931.68 



Jarque-Bera Test: We implemented the Jarque-Bera test to check the null hypothesis Ho 
that the clearing prices ir(t) are normally distributed (while Hi is the hypothesis that n(t) are not 
normally distributed). The Jarque-Bera test used both skewness and kurtosis simultaneously to 
check the normality of the selected data. The Jarque-Bera test statistic was calculated as 

JB = (S*) 2 + (K*) 2 , 

where 

s * = \/i § {7r{t)) ~ iV(0 ' 1} and K * = \f^(* {7T{t)) " 3 ) ~ iY(0, 1} - 

The Jarque-Bera test yielded the p- value of 3.161%, which is below the default significance 
level of 5%. Thus, we rejected the null hypothesis Hq that the clearing prices ir(t) are normally 
distributed. However, with the p- value of 3.161%, the distribution of n(t) is rather close to a normal 
distribution, which is also shown in Figure 1. 




\ \ \ 1 1 1 1 

20.0 20.1 20.2 20.3 20.4 20.5 20.6 
Clearing prices (US dollars) 



Figure 1: The density of clearing prices n(t) 
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Evolution of ir(t): From the result of linear 
drift of clearing prices vr(t). Thererfore, we have 

dir(t) = cdt + cr w (t) / b w 

J s 



regression, we obtained the constant value of the 

(s,t)W(ds,dt), 
where: c = 0.000374764. 




Buy Limit Orders: Let us recall the definition of Q 

log Q(—K, iAt) = log number of buy orders * buy order quantity arriving before iAt^J , 

where At = 1 minute. 

We tested the data of logQ(iAt) for one-minute time interval with the autocorrelation and 
partial autocorrelation functions. The ACF and PACF figures showed (see Figure 3) that the data 
of log Q (iAt) for one-minute time interval were fitted in the autoregressive AR(1) model. 





Lao 



Figure 3: ACF and PACF of log Q (iAt) for one-minute time interval. 



Results: We conducted the estimation of our model (18), (19) over the data set described above. 
The opening clearing price on that day was 

?r(0) = 20.16. 
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Using K = 7 and Ap = 0.05, the speed of mean-reversion was obtained from the autoregressive 
AR(1) model. For the relative net demand curve, we observed the following initial value and relative 
hourly volatility: 



Q(-K,0) = 1.02705 x 10 11 , 
(Jq\—K) = 0.01976. 

For the relative net demand curve density the values are reported in the following table 



k 


q{k,0) [in 10 11 ] 


a rei (k) - hourly 


a{k) 


-6 


0.95314 


0.04883 


0.11903 


-5 


1.41994 


0.04655 


0.29142 


-4 


2.35893 


0.01706 


0.25250 


-3 


0.82541 


0.04423 


0.36708 


-2 


0.14050 


0.04877 


0.36752 


-1 


4.13487 


0.03744 


0.29380 





0.21397 


0.00461 


0.21991 


1 


9.95599 


0.00379 


0.25219 


2 


4.61052 


0.03496 


0.36316 


3 


3.51037 


0.00653 


0.15248 


4 


2.42507 


0.01224 


0.19830 


5 


0.14219 


0.00036 


0.34405 


6 


2.70257 


0.00969 


0.13387 



Table 2: Statistics for the relative net demand curve density 



We then simulated the lognormal model (18), (19) using N = 100 scenarios and approximated 
the value of a call option: 



C ^ = jr E max (^ ") - v>, o), 

where ip is a strike price for an expiration T=0.02, that is, roughly equal to one week. We then 
calculated the implied volatility of the call for each strike price, that is, the value a tmp (Strike) such 
that the Black-Scholes value of the call option equals C(K), and this for each K. The resulting 
function a imp (Strike) as a function of the strike price ( the smile curve) is reported in Figure [4] 
below. 

As expected, the smile is fairly pronounced, which is an indicator that our model engenders fat 
tails in the risk- neutral distribution of 7r(T). 
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Figure 4: Implied volatility as a function of the strike price. Call option with 1 week expiration, 
7r(0) = 20.16, and zero interest rate. 



5 Appendix 

Lemma 5.1. Let 

fip(x,t) = 

ap(x,t) = 
b P (x,s,t) = 



V Qa (Pa(x, t), t) + l^(P A ( x , t),t)a 2 P (x, t) + ^(Pa(x, t),t)a p (t) 



§(P A (x,t),t) 



VQ A ( P A(x,t),t) 
^(P A (x,t),t) 
bQ A (PA(x,t),t),s,t) 



Then P A is a semimartingale and satisfies 

S 

dPA(x,t) = fj,p(x,t)dt + ap(x,t) j bp(x, s,t)W(ds,dt) 



Proof. By definition: 

Q A (P A (x,t),t) = x 



(20) 



(21) 



We suppose that (20) holds and apply the Ito-Wentzell formula (see e.g., |Kry 09l) to both sides 
of pi} yields: 



^Q A ( P A{x,t),t) + -Q—(P A (x,t),t)fl P (x,t) 

+ l^(P A (x,t),t)a 2 P (x,t) + d ^(P A (x,t),t)a P (x,t) = 0, 
VQ A (PA{x,t),t)bQ A (P A {x,t),s,t) + -^-(PA(x,t),t)a P (x,t)bp(x,s,t) = 



□ 



Proof of Theorem 2.2 
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For fact (Fl), we first check assumptions (RF) in |KR09j : 

(i) Since Lz,(i?, t) is a semimartingale, it is a strong integrator. 

(ii) We calculate 

^,t)_a _ i 



a» 2 a* "" %(p A ({,.t)) 

Since Q A is strictly decreasing and twice differentiable in p, then its inverse Pa is differentiable 
in i9 and - L gJ^'^ is continuous in i?. 

(iii) Since ^-^sr^ = PA($,t) is a semimartingale, it is a strong integrator. 



iv) By the Lemma 5.1 the quadratic variation of ■) is equal to 



This is clearly strictly increasing if (tq(P a (x, t), t) is uniformly bounded away from zero, which 
holds if ctq(p, t) is uniformly bounded away from zero on < p < S. 

Thus assumptions (RF) are satisfied. We now define the market price of risk as 

A (. )( ) = MM) 

Since p,p(x,t) is bounded and ap(x,t) is non-zero, then the assumption (UB) in |KR09j is 
satisfied. To check the assumption (UI) we refer the reader to |KR09| for the definition of A* and 
Since both Qa and Ql are uniformly bounded, the integrand is bounded. By the 
same reason, and because of the continuity (in time) of Qa, the quadratic variation [M^ n \ M^] is 
bounded. Thus (UI) holds, and Theorem 3.5 in [KR09] holds, proving Fact (F2). 

Fact (F2) follows directly from the comments after Theorem 3.5 in [KR09J, which we quote here, 
while adjusting for our notation and formulae numbering: 

"Consider now again the dynamics ^ of the real wealth process, and let 6 be such that 
f Ll(9,(1s) is bounded from below (the transaction costs term in ^ can be avoided by the large 
trader by using only tame strategies). An arbitrage opportunity is an admissible strategy such 
that we have for the associated real wealth process V e that V e (0) < 0, V s (T) > P-a.s., and 
P(V (T) > 0) > 0. By Theorem 3.5 in |KR09j . there exists a probability measure Q such that 
f Ll(8, ds) is a Q-local martingale, hence a supermartingale. It follows now from the dynamics ([5]) 
of the real wealth process that Eq[V (T)] < V e (0) which, as Q is equivalent to P, excludes arbitrage 
opportunities for the large trader. " 

We can now prove the continuity of vr(t). Since by assumption uncross orders result in Ql being 
continuous in time, a discontinuity can arise only from cross orders. Suppose that the large trader 
is a net buyer (for a sell strategy the argument is identical), i.e, that there is an e > such that for 
all < 5 < e 

9(t + 5)-6{t) > 0. 

Suppose that 

7r(t) -7r(t_) > e 2 > 0. (22) 
Then, since 8Qa/9p > and Qa is assumed to be twice differentiable 

-Q A (*(t),t) + Q A (^(t-),t) > e 3 > 0. 
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By continuity in time of Q A 



-Q A (n(t),t + ) + Q A (Tr(t-),t) > e 4 > 0. 



An examination of terms in Lemma 3.3 in ! >!>() 1 shows that discontinuous strategies are suboptimal 
for the trader, i.e. 



By definition 



0(t+) - 9(t) = 0. 
9(t+) - 6(t) = -Q A (n(t),t + ) + Q A (ir(t-),t). 



(23) 



However, this shows, that ( |22[ ) contradicts (23). Thus ir is continuous, proving fact (F4), and 

Q L (n(t),t)-Q L (ir(t-),t) = 0. (24) 

By definition also 

0(t+) - 6(t) = Q L (Tr(t),t + ) - Q L (*(t-),t) (25) 
= Q L (7r(i), t+) - Q L ^{t),t) + Q L (ir(t),t) - Q L {*{tJ), t). 



Combining (23), (24) and (25) we obtain that 

Q L (ir(t),t+)-Q L (Tr(t),t) = 0. 

This shows that the net demand curve is continuous in t, thus proving fact (F3). 
Finally, we build the curve: 

Q = Ql + Qa- 

Suppose that a new atomistic trader comes to the market. She will trade at a price 7r(i) and 
will have liquidity costs: 

t 

rt 



j L(0(u),du) = J 9(u)dir(i 



Since J L(6(u),du) is a Q-martingale, therefore it is a Q-martingale, too. 
o 



□ 
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